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Abstract 

I introduce a new approach to the maximum flow problem by a 
simple algorithm with a slightly better runtime. This approach is 
based on sorting arcs insight of vertices on a residual graph. This 
new approach leads to an O(mn 5 ) time bound for a network with n 
vertices and m arcs. 

The Category: Algorithms, Graph Theory and maximum flows 

1 Introduction 

The maximum flow problem and the minimum cut problem have a big variety 
of uses in scientific and engineering applications. Even if it's one of the well 
understood combinatorial problems with a long history and many approaches 
to improve the asymptotic runtime. 

In this paper an algorithm which combines calculation of blocking flows, 
shortest path and the usage of preflow[2], will be introduced. By intensively 
studying the hi _prf the author wanted to develop a non numerical topology 
based algorithm to get rid of the relabel function. Coming up with the 
idea, that it might be possible to know where the excess needs to move, can 
be calculated faster. For now that can be done, but not for a single vertex. 
More for using something similar to a global update [3 J to calculate something 
similar to an acyclic core, known from Binary Blocking Flow pQ. 
The advantage is the speed that core can be calculated in 0(m), even though 
that calculation is very simple it generates all necessary informations, to know 



all dependencies for each residual network, for minimizing the amount of non 
saturating pushes. Though there is only a global update, which is using unit 
distance. 

During the global push all improving dependencies will be augmented. As 
shown later this will also only need 0(m), because every arc of the core will 
be only once augmented. 

Even with worst case improvement of each iteration it can be shown, that 
those two functions only need to be called 0(n°' 5 ). 

2 Definitions and Notation 

The input to the maximum flow problem is a network N(G(V, E), s,t,u), 
where s, t e V : s ^ t are source and sink and u : A — > N is the capacity 
function. Let n = \V\ and m = \A\. 

A flow in a network is a function / : E — > N where for each arc, the capacity 
constraints < f(a) < u(a) holds and for each vertex j e V / {s,t} the 
conservation constraint Y,(j,k) f{h — f(hj) = holds. The value of 

the flow is \ f\ = J2 (j,t)f(j,t). f(a) is referred as the flow on arc a. 
A preflow is relaxing the conservation constraint to be J2(j,k) f{j-> k) < Y^nj) f(hj) 
resulting in an excess e(j) = Yhj) f(hj) — J2(j,k) /C?\ k). We assume that 
every arc a has a reverse arc a R . The residual capacity Uf of an arc is 
defined to be u(i,j) - f(i,j) + f(j,i). 

Let o be a function \/v3o : N — > (v , w) E A order of outgoing arcs for each 
vertex. 

3 Sorting Flow 

The algorithm has only two operations effecting the whole graph. The first we 
call Breadth First Search Sort or short BFSS, which calculates the arcs which 
need to be augmented. The second is a global push operation augmenting 
these arcs. These two operations are repeated until no augmenting path can 
be found. BFSS will return true, when an augmenting path exists. 
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Require: N(G(V, E), s, t, u) 
Ensure: |/| to be maximized 
1: while BFSS do 

2: Push() 
3: end while 

Algorithm 1: Sorting Flow 

3.1 BFSS 

This function has two objectives. First: checking for an augmenting path, 
which here means vt : e(v) > 0, where Va e vt : Uf(a) > 0. The second 
objective is to sort all arcs in every vertex by their shortest distance. To 
achieve this we run a regular BFS starting with t, using a queue and looking 
at ingoing arcs. During this every vertex has one of three states: not found, 
1 found and 2 used. Every vertex in state 1 is element of the queue. 
In initial all nodes are state 0, except t, which also is in the queue. While the 
que isn't empty the algorithm removes the first element of the queue, marks 
it as state 2 and now iterates through the ingoing arcs (u, v). For all vertices 
u which state is 0, set their state to 1 and queue them to the queue. For all 
vertices with a state ^ 2 set the next free value in o(u) to the arc (u,v). 
When ever a vertex v is removed from that queue, if e(v) > or v = s, put 
this node into another queue, to not mess up those to queues the second one 
will be referred as priority queue. 

At the end of this function we have a priority queue, with all vertices e/(i>) > 
0. Now looking only on arcs addressed by function o reveals an acyclic core. 
The ordering insight of the priority queue matches for every residual graph 
with with acyclic core. 
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1: final = true; 

2: que.push(t) 

3: while !que.empty() do 



4: v — que. get Front () 

5: if e(v) > then 

6: priority queue. pushLast(v) 

7: end if 

8: f.state=2 

9: for all u f (u,v) > G E do 
10: add (u, v) to o(u) 

11: if u. state 7^ 2 then 

12: que.pushLast(w) 
13: -u.state=l 
14: end if 

15: end for 



16: end while 

Algorithm 2: Breadth First Search Sort 



3.2 Push 

In this step the excess will be pushed as much as possible. By using the 
priority queue and only augmenting arcs, which are addressed by the function 
o. As long as the priority queue isn't empty the algorithm will remove the 
first node v and will as long as there are arcs in o(v) and excess e(v) > 0. All 
arcs of o(v) will be augmented. If the the excess of another vertices changes 
from to > 0, those vertices will be put to the front of the priority queue to 
be executed next. 



1 


while !priority queue. empty () do 


2 


v =priority queue.getFront(); 


3 


for all (v,w) G o(v) or until e^iy) = do 


4 


if ef(w) = then 


5 


priority queue. pushFront(w) 


6 


end if 


7 


augment (v, w) 


8 


end for 


9 


end while 



Algorithm 3: Push 
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4 Correctness and Runtime 



4.1 Correctness 

Even in a worst case situation BFSS will find at least one augmenting path, 
if this exists. This path will be found in Vo(w) £ vt : e.f(y) > 0. The Push 
will increase the excess of the sink at least by one and saturate one arc or 
discharge one vertex. If no augmenting path vt : ef(v) > can be found the 
algorithm determines. According to Ford and Fulkerson this blocking flow is 
supposed to be maximal. 

4.2 Runtime 

BFSS and Push look at every arc only once, therefor the runtime for each 
call of these functions has to be 0(m). 

The amount of calls for those two functions can be easily understood, when 
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Figure 1: to be solved in one iteration 

thinking about the smallest worst case inputs. Starting with the graphs, 
which can be solved in one iteration. Those can be described as a simple line 
as seen in figure [1} To force a second iteration it is necessary to block the 
flow at one arc, as seen in figure [2j If d > c > b > a the excess will stop on 
the middle layer of this graph. After the next iteration is done. The max 
flow will be calculated. 

This can be done for even one layer more or an unlimited amount of layer as 
seen in figure [3j When deforming these graphs, and observing the algorithm 
in their worst case it still can be observed, that the algorithm completes one 
layer breadth and one layer length in each iteration, concluding to an amount 
of calls 0(min(length(f),breadth(f)) C 0(n°- 5 ), where / is the calculated 
maximum blocking flow. 
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Figure 2: to be solved in two iterations 




Figure 3: to be solved in three and n — 1 iterations 
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